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Abstract 



We introduce the matrix sums that represent a discrete analog of the 
matrix integrals of random matrix theory. The summation runs over the 
set T n of all possible n-vertex graphs 7„ weighted by exp{— f3 TrA n }, (5 > 
0, where A„ = A(7„) is the analog of the Laplace operator determined on 
j n - Corresponding probability niCfiSLirG on L n reproduces the well-known 
Erdos-Renyi ensemble of random graphs. Here it plays the same role as 
that played by the Gaussian Unitary Invariant Ensemble (GUE) in matrix 
models. 

Regarding an analog of the matrix models with quartic potential, we 
study the cumulant expansion of related partition functions. We develop 
a diagram technique and describe the combinatorial structure of the co- 
efficients of this expansion in two different asymptotic regimes ft — O(l) 
and /3 = 0(log n) as n — » oo. 

Keywords: Laplace operator on graph, partition function, Erdos-Renyi random 
graphs, connected diagrams, Catalan numbers, Polya equation. 

1 Introduction 

During last three decades, the studies of matrix models of theoretical physics 
have deeply influenced a number of branches of modern mathematics and math- 
ematical physics. The central notion here is the partition function Zn((3,Q) 
given by the integral over the set Mn of all Hermitian iV-dimensional matrices 



where Q is a "potential" function, dH is the Lebesgue measure, Cjy = Zn((3, 0) 
is the normalizing constant, and E^ E {-} denotes the mathematical expectation 
with respect to the probability measure Vn with the density 



H 




cxp{-/3Trff 2 + Q(H)}dH = C N E^ E {e Q{H) }, (1.1) 



N 



exp{-/m-ff 2 } = exp{~/3 £ \H V f}. 



(1.2) 



The measure Vn supported on A4n generates the Gaussian Unitary Invariant 
Ensemble of random matrices abbreviated by GUE (see monograph [12] for the 
detailed description of the ensemble and its properties). 

The first non-trivial example of (1.1) is given by the quartic potential 

Q(H) = jjTrH 4 . (1.3) 

The matrix model (1.1), (1.3) has served as the source of a series of deep results 
establishing connections between orthogonal polynomials, integrable systems, 
moduli spaces of curves and such combinatorial structures as maps (see [3] and 
[5] for the earlier and more recent results and references and [4] for the review). 
One of the simplest but important result is that the leading term of the formal 
asymptotic expansion of variable 

^logE^jexp^Trtf 4 )} (1.4) 

in the limit N — > oo is given by the series in powers of g with the coefficients 
determined by the numbers of 4-valent two-vertex maps dual to the famous 
quadrangulations of the compact Ricmann manifold (see [19] for introductory 
description and references therein). 

It should be noted that the GUE and its real symmetric and symplectic 
analogs represent very special class of random matrices. It is natural to ask 
about matrix models when the mathematical expectation in (1.4) is taken with 
respect to a measure different from that determined by GUE. 

In present paper we introduce a discrete analog of (1.1), where the integra- 
tion over at is replaced by the sum over the set T n of all possible simple graphs 
7„ with n vertices. In this setting the "kinetic energy" term TrH 2 is replaced by 
TrA„, where A„ = A(7„) is the discrete analog of the Laplace operator deter- 
mined on the graph 7„. Corresponding Gibbs weight exp{— TrA(7„)},/3 > 
generates the probability measure n n {0) on T n . 

A simple but non-trivial property of this ensemble plays a very important 
role in what follows. The observation is that the probability space (T n , /j, n ) 
coincides with the widely known Erdos-Renyi ensemble of random graphs with 
n vertices (see e.g. [2]). In this ensemble, the indicator functions of edges 
are represented by jointly independent Bernoulli random variables. As far as 
we know, this connection between Erdos-Renyi random graphs and the Gibbs 
measure C ra exp{— TrA(7„)} was not observed before. 

Our aim is to consider the asymptotic behavior of F n = logE Mri {exp(Q„)} 
where Q n is given by the analog of quartic potential (1.3) and to explore the 
combinatorial structures that arise in this problem. We develop a diagram 
technique to study the cumulant expansion of F n . We show that the leading 
terms of this expansion are related with the number of connected diagrams on 
the set of two-valent vertices. We derive recurrent relations for the numbers 
of such diagrams and describe the coefficients of the cumulant expansions of 
F n in two different asymptotic regimes when = 0(1) and — O(logra) as 
n — > oo. These recurrent relations generalise those for the Catalan numbers. 
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Corresponding generating function verifies an equation similar to the Polya 
equation for the generating function of the rooted Cayley trees. 



2 Graph Laplacian, matrix sums and random 
graphs 

Given a finite graph with the set V n = {v\, . . . ,v n } of labelled vertices and 
E m = {e {1 \ . . . , e(™)} the set of simple non-oriented edges, the discrete analog 
of the Laplace operator A(7) on graph 7 can be determined as (see for example, 
[14]) by relation 

A( 7 ) = d*d, (2.1) 

where d is the difference operator determined on the space of complex functions 
on vertices V n — > C and d* is its conjugate determined on the space of complex 
functions on edges E m — > C. 

It can be easily shown that in the canonical basis, the linear operator A(7) = 
A n has n x n matrix with the elements 

(deg(vj), if i = j, 
Ay = < -1, if i¥= 3 and (v it Vj ) G E, (2.2) 

[ 0, otherwise, 

where deg(w) is the vertex degree. If one considers the n x n adjacency matrix 
A = ^(7) of the graph 7, 



A = ( 1, if {vuVj) G E,i^ j, 
■ 0, otherwise, 



(2.3) 



then one can rewrite the definition of A (2.2) in the form 

n 

Ay = Bij - Aij, where B l3 = 5 tj ^ A iU (2-4) 

1=1 

where 5 is the Kronecker symbol 

r _ / 1, if i = 3, 

" I if^i- 

It follows from (2.1) that A(7„) has positive eigenvalues. 

Let us consider the set T n of all possible simple non-oriented graphs 7„ with 
the set V = V n of n labelled vertices. Obviously, |r„| = 2"(™- 1 )/ 2 . Given 
an element 7 e r„, it is natural to consider the trace TrA(7) as the total 
energy of the graph 7. Then we can assign to each graph 7„ the Gibbs weight 
exp{— (3 TrA(7„)}, (3 > and introduce the partition function 

Z n (P,Q)= exp{-/3 TrA„ + Q( 7 „)}, (2.5) 
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where A„ = A(7„) and Q is an application r„ — » R that we specify later. Using 
the fact that TrA = Trd*d, one can consider (2.5) as a discrete analog of the 
partition function (1.1). We will see that this is especially interesting in the 
case of quartic potential (1.3). 

Let us note that we should normalize the sum (2.5) by |r„|, but this does 
not play any role with respect to our results. In what follows, we omit subscript 
n in A„. 

Definition (2.4) implies that 

n n 

TrA = ^A ii= £ ^=2 E ( 2 ' 6 ) 
Then we can rewrite (2.5) in the form 

7n£r„ l<i<j<n 

It is easy to see that 

Z n (P,0) = (l + e-^) n(n - 1)/2 . (2.8) 
Then the normalized partition function can be represented as 

Z n ((3,Q) = Z n (f3,Q)/Z n ((3,0) = E p {e^} , (2.9) 

where E^j-} denotes the mathematical expectation with respect to the measure 
supported on the set T n . This measure assigns to each element 7 G T n the 
probability 



^(7) 



e -2/3|S(7)l 



^ + p -2/J)™(™- 1 )/ 2 : 



where E{~j) denotes the set of edges of the graph 7. 

Given a couple (x, y), x,y € {1, . . . , n}, one can determine a random variable 
a xy on the probability space (T n ,P n ) that is the indicator function of the edge 

{Vx,Vy) 

1, if (v x ,v y ) e E{i), 
0, otherwise. 



It is easy to show that the random variables {a xyi 1 < x < y < n} are jointly 
independent and are of the same Bernoulli distribution depending on (3 such 
that 

a (« = J X ' with Probability ^—jp = P, (2 10) 

t 0, with probability 1 — p. 

The probability space (T n ,P n ) is known as the Erdos-Rcnyi (or Bernoulli) 
ensemble of random graphs with the edge probability p. Since the series of 
pioneering papers by Erdos and Renyi, the asymptotic properties of graphs 
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(r„, P n ), such as the size and the number of connected components, the maximal 
and minimal vertex degree and many others, are extensively studied (see [2, 7]). 
Spectral properties of corresponding random matrices A (2.3) and A (2.4) are 
considered in a scries of papers (in particular, see [1, 6, 8, 9, 10, 13, 16]). In 
present paper we study the random graph ensemble (r„, P n ) from another point 
of view motivated by the asymptotic behavior of partition functions (2.9) with 
certain "potentials" Q. 

3 Partition functions and diagram technique 

In the previous section, we have shown that the Gibbs weight exp{— [3 TrA(7)} 
generates the probability measure on graphs equivalent to that determined by 
the Erdos-Renyi ensemble of random graphs. This gives us an important tool 
for direct computation of averages of the form (2.9). Another important point 
is that this Gibbs weight leads us to the correct definition of discrete analogs 
of the matrix models (1.1) and in particular of the matrix models with quartic 
potentials (1.4). That is why one cannot neglect the Laplacian form of the 
Gibbs measure and start simply with computations of averages with respect to 

Pn- 

3.1 Analog of the quartic potential 

Let us determine the discrete analog of the partition function (1.1) with quartic 
potential (1.3). Once Tr_ff 2 replaced by Tr(9*<9) = TrA, it is natural to consider 

Tr(d*dd*d) = TrA 2 

as the analog of TrH 4 . Then the partition function (2.5) reads as 

Z n {p,g)= ]T cxp{-/?TrA„+.g„TrA 2 }, (3.1) 

7^er„ 

where g n is to be specified. It follows from (2.3) and (2.4) that 

n n 

TrA 2 = TiB 2 + TrA 2 = ^ (A 2 ) t] + ^ A l0 . 

Then, using (2.6) and repeating computations of (2.7) and (2.8), we obtain 
representation 

Z n (0,g) = Z n (f3,g)/Z n (f3,0) = f ^—^ J E /3 /{e 9 " x "}, (3.2) 

In this relation, we have denoted (3' = (3—g n and introduced the random variable 

n 

X n = ^ auaij, (3.3) 
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where are jointly independent random variables of the law (2.10) with 
replaced by 1 '. The average E^/ denotes the corresponding mathematical ex- 
pectation. In what follows, we omit the subscripts and 0' when they are 
not necessary. We study the limiting behavior of the cumulants of the random 
variable X n in two different asymptotic regimes. The first limiting transition 
is determined by the choice n = O(logn) as n — > oo. We study this case in 
Section 4. In Section 5, we consider the second asymptotic regime given by of 
n — const. 

3.2 Cumulant expansion 

Using the fact that Q>ij aPC bounded by 1 and X n < n 3 , we can write that 

OC 

logE{e»*»} = J2^Cum k (X n ), (3.4) 

k=l 

where Cumk(X n ) is the k-th cumulant of random variable X n . 

Cum k (X n ) = A_ (l og E{e^"}) | fl=0 . (3.5) 

Denoting by Y a the random variable auaij with triplet a = (i, I, j), we can write 
that _ 

Cum k (X n )= ^2 cum{Y ai , . . . ,Y ak }, (3.6) 

{ai,...,a k } 

where the sum runs over all possible values of a s , s = 1, . . . , k and 

d k 

cum{Y ai , . . . , Y ak } = — logE{exp(zir Ql + . . . + z k Y ak )}\ Zr=0 . (3.7) 

azi ■ ■ ■ az k 

The variable cum{Y ai , . . . , Y ak } is also known in probability theory as the semi- 
invariant of k random variables Y. 

Let us introduce a graphical representation of the set of variables Y. Given 
k values ol\, . . . , a k , we represent the set of random variables Y ai , . . . , Y ak by 
the set of k labelled vertices U = (m, . . . , u k ). Regarding two variables Y Ur and 
Ya s , we join corresponding vertices u r and u s by an edge s(r,s) if and only if 
Y ar and Y ag have at least one variable a in common. Considering all possible 
couples (r,s),l < r < s < k and drawing corresponding edges, we obtain a 
graph that we denote by G k = (U k ,£ k ). It is clear that this graph G k = G k (a k ) 
depends on particular value of the variable a k — {ct\, . . . , a k ). The following 
proposition is a well-known fact from the probability theory of random fields 
[11]. 

Lemma 3.1. The semiinvariant cum{Y ail . . . , Y ak } is not equal to zero if and 
only if the graph G k (a k ) is connected. 

Proof. Let us first note that if two random variables Y Ur and Y" Qs have no 
variables a in common, then they are independent. Clearly, the same obser- 
vation is true for the subfamilies of random variables Y. This means that if 
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the graph Gfc(dfc) consists of two or more non-connected components, then the 
corresponding subsets of random variables Y are jointly independent. The char- 
acteristic property of the semiinvariants is that it vanishes in this case [11]. This 
completes the proof of Lemma 3.1. 

Let us note that the fact that Y ar and Y ag have one or more variables a in 
common means that corresponding variables i,l,j coincide. In particular, the 
fact that Y ar and Y as have exactly one variable a in common implies that one of 
the eight possibilities for the sets (i r ,l r ,jr) and (i s ,l s ,j s ) occurs. For example, 
this happens when 

[i T is^r Is? jr js) or {i r l s , l r i s , j r ^ js)- (^-^) 

We will refer to such cases as to the direct and inverse gluing, respectively. 

Lemma 3.2. Given /3 and g fixed, the number of terms in Cum,k{X n ) is given 
by relation 

#{Cum k (X n )} = 0{n k+2 ) asn^oo. (3.9) 

Proof. To draw an edge e(r, s) of the graph Gk means to make equal at least 
two variables a taken from Y ar and Y as , respectively. To make equal two or 
more a's means to make equal some of the variables (i r , l r ,j r ) and (i s , l s ,j s ). 

Let us describe the process of drawing edges of G/. step by step. One starts 
with 3fc variables i, j, I that can take values 1, . . . , n independently. Each gluing 
of two variables a diminishes by 2 the number of variables that move indepen- 
dently. To make the graph connected, we have to draw at least k — 1 edges, so 
we are forced to perform at least k — 1 gluings. When this is done, the number 
of variables that move independently is less or equal to 3k — 2(k — 1) = k + 2. 
Maximizing the number of variables that can take different values, from 1 to n 
independently, we obtain 

n(n — 1) ■ ■ ■ (n — (fe — 1) + 1) = 0(n k+2 ) 

terms. Lemma 3.2 is proved. 

3.3 Connected diagrams 

Let us further develop the graphical representation of the set of variables Y and 
give more details for the description of corresponding connected graphs Gk ■ At 
this stage we make no difference between the direct and inverse gluings (3.8). 

First let us note that each variable Y a by itself can be represented by a 
graph of three vertices corresponding to variables i, I and j joined by two edges 
denoting random variables an and aij , respectively. Slightly modifying this, we 
can say that Y ar is represented by a vertex u r with two off-spreads representing 
variables a ir i r to the left and ai r j r to the right from u r . We say that the vertices 
u r are two-valent. In Gk, there are k such two-valent vertices joined by k — 1 
edges that we refer now to as to arcs. The arcs join two different off-spreads 
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that wc glue during the procedure described in the proof of the Lemma 3.2. 
The family of k two-valent vertices together with k — 1 arcs represent a diagram 
that we denote by 5k- This diagram provides more information than the graph 
Gk because it shows exactly which off-spreads are glued between themselves. If 
one forgets the two-valent structure of u r 's, one gets the tree Gk- One of the 
possible example of 8k and corresponding Gk for fc = 5 is given on figure 1. 




1 2 3 4 5 12345 



Fi gure 1: Connected diagram 6^ and corresponding tree G k for k = 5. 

One can see that several off-spreads can be glued together. In this case 
we say that they are colored by the same color. We call the off-spreads that 
remain non-glued as the free off-spreads and leave them non-colored (or grey). 
On figure 1 the diagram 5$ contains three grey elements and three color groups 
of 2, 2, and 3 elements. 

Clearly, one can draw several diagrams that represent the same coloring of 
the off-spreads. An example of such two diagrams 5k and 5' k with k — 5 is given 
on figure 2. Wc call the diagram the reduced diagram when the off-spreads of the 
same color are joined by the arcs connecting the nearest neighbors. On figure 2 
the diagram 5k is the reduced one, the diagram 5' k is not. 




1 23451 2345 



Figure 2: Non-reduced diagram S' k and the equivalent reduced one. 

This argument explains the difference between the diagrams we have and the 
set of trees on labelled vertices. Indeed, the diagrams 5k and 5' k are equivalent 
but the corresponding trees are not. 

With this diagram representation, we see that the leading contribution to 
(3.9) comes from the family T>k of connected acyclic reduced diagrams 5k drawn 
on the set of k two-valent vertices with ordered off-spreads. The number of such 
diagrams dk — \T>k\ is determined in the next section. We complete this section 
with the following simple proposition. 

Lemma 3.3. Let us consider a connected reduced diagram 5k that have r 
color groups of arcs, with fj,i, jj,2, ■ ■ ■ , /U r arcs in each color group. Then 

Mi + M2 + • • • + Mr = k - 1. (3.10) 
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There are k — r + 1 grey elements in 5k ■ 

Proof. Each diagram generates a tree on k vertices u\, . . . , Uk- So the total 
number of arcs is equal to k — 1. The group of /i s arcs of the same color produces 
fi s + 1 color elements. The total number of colored elements is k — 1 + r. Then 
the number of non-colored (grey) elements in 5k is 2k — (k — 1 + r) = k — r + 1. 



4 Sparse random graphs 

In present section we study the case when the "temperature " T n = l//3„ van- 
ishes when n — > oo 

Pn = \ log Q) (l + o(l)), n^oo (4.1a) 

with some c > 0. This corresponds to the random graph ensemble (3.2) with 
vanishing edge probability 

Pn = y-J , n -> oo, (4.16) 

where we denoted c = ce 2g ^ n ^ and g(n, c) = 5/c. In this paper we consider the 
limiting transition when n, c — > 00 and c = o(n). The case of n — > 00, c = const 
will be studied in separate publication. The main results of this section are 
presented by the following two statements. 

Theorem 4.1. 

Given k > 2, there exists the limit 

lim -L-Cum k (X n ) - 2 fe - 1 d fc , (4.2) 

ra,c->oo,c=o(n) UC k+L 

where the numbers dk , k > 2 are determined by following recurrent relation 

d k = 2kd k -i + ( • ) (i + l)(fc - i) djdk-i-j , k>3 (4.3) 

with the initial condition d 2 = 4. 

Let us consider the auxiliary numbers 

h k = ^± for k>2. 
k\ 

It is easy to deduce from (4.3) that the sequence h can be determined by the 
following recurrent relation 

k + l ^ 

h k = 7 2_/ ^ hk-i-j, h a = 1. (4.4) 
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Proposition 4.2. 

The generating function h(z) = X^fcLo hkZ k is determined in the complex domain 
Ri/8 = {z € C : \z\ < 1/8} and verifies there the following equation 

h(z) = cxp{2z/i(z)}. (4.5) 

It follows from (4-5) that 

h k = 2 k ik + V k ~\ k>l. (4.6) 

Using (4.4), it is easy to find the first values of d 

d 2 =4, d 3 = 32, rf 4 = 400, d 5 = 6912, d 6 = 153664. 

The general form of {d k } follows from (4.6) 

d k = 2 fe (fc + l) fe - 2 , k>2. (4.7) 

Equation (4.5) is similar to the Polya equation for the generating function of 
the rooted trees on labelled vertices [15, 17] (see also relation (4.17) below). The 
explicit expression (4.6) resembles the number of Cayley trees on k + 1 vertices. 
However, we did not find any obvious one-to-one correspondence between the 
set of Cayley trees on k + 1 vertices and the fc-vertex diagrams we count. While 
the skeletons of our diagrams are given by trees, the maximal degree of these 
trees is bounded by 4. Also, there are equivalent diagrams that have different 
tree skeletons. This makes the set of the diagrams we study quite different from 
the family of Cayley trees. Relation (4.4) generalizes recurrent relations for the 
Catalan numbers. Up to our knowledge, the class of diagrams V k as well as the 
numbers d k were not considered previously. 

Proof of Theorem 4.1. Let us separate the set of all possible values of triplets 
oji,. .. ,a k containing n 3k elements into the classes of equivalence labelled by 
5h € T>k- This is done in obvious way. Then we can write that 

cum{Y ai , . . . ,Y ak } = Y Af(S k ) cum{Y(8 k )}(l + o(l)), n -> oo 

ai,...,a k S k eV k 

(4.8) 

where M{8 k ) denotes the number of elements in the equivalence class labelled 
by S k , 

cum{Y(5 k )} = cum{Y a(Sk)l Y a(Sk)k }, 

and a(S k ) is one of the representative of this equivalence class. For instant, we 
can choose a(5 k ) = (a\, . . . ,a k ) with minimal possible values of i\, l\, ji, 12, 
and so on, in the way such that (ai, . . . , a k ) belongs to the class S k . It follows 
from Lemma 3.2 that 

M{5 k ) = n{n -!)••• (n - k- 1) = n k+2 {\ + o(l)), n 00. (4.9) 
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Using the basic property of the semi-invariants, we can write that [11] 

cum{Y(5 k )} = £ E{F Tl( ^ (5fc)) } • • • E{Y TA7rkMSk)) }(-iy-\a-l)\ , (4.10) 
ft 

where the sum runs over all unordered partitions T\ , . . . , T„ of the set of vertices 
Uk = {u\, . . . , Uk}; that is over all families of non-empty non-intersecting subsets 
T s of Uk giving in sum all Uk- The second argument of T's reminds us that the 
variables a(5k) s attached to u s are chosen according to the rule prescribed by 8k- 
Random variable Y T r nk . iBi /g k \\ is given by the product of corresponding random 
variables Y. 

Let us consider the term of the sum (4.10) that corresponds to the trivial 
partition 71^ of Uk consisting of one subset only: 7\ = Uk- This term is given 
by the average value 

E{F^ )l ...y s(5fe) j = w(4 1) ). 

According to Lemma 3.3, there are k — r + 1 non-colored elements in diagram 
5k, where r is the number of color groups of arcs in 5k- These grey elements rep- 
resent independent random variables o that are also independent from elements 
belonging to color groups. Then we obtain the factor 

(Ea) k - r+1 = (-) 

Each color group produces the factor Ea^^ 1 = Ea = c/n and there are r color 
groups. Then the weight of the partition ir^ is 



^(4 1) ) = (r) • (4- 11 ) 



fe+i 



Now it is clear that any other partition 7Tfe produces the term of the order 
o((c/n) fe+1 ) that is evidently smaller than that of (4.10). More precisely, we can 
write that 

/ c\ k+1+x 

E{Y Tl(nk . MSk)) }---E{Y TA n kM s k ))} = [-) , (4.12) 

where \ lSi equal to the number of the arcs of 5k such that the left and 
right feet of 9 belong to different subsets T' and T" of the partition Uk under 
consideration. We say that these arcs are cut by partition 7^. In particular, the 
partition tt^ of Uk into fc subsets produces the factor 



■*<•?'> -(h) 



2A- 



because all k — 1 arcs of 5k are cut by this partition. 

Remembering that each arc can be drawn in the direct and inverse sense 
(3.8) and gathering formulas (3.6), (4.8), (4.9) and (4.11), we conclude that 

-4— -Cum k (X n ) = 2 fe - 1 4(l + o(l)), l«c«n, 
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where dk = \T>k\ is the total number of diagrams 6k- Let us derive recurrent 
relations for dk- 

Let us consider one particular diagram 6k and denote by Pj(6k) some part 
of 6k that consists of certain j vertices and all arcs that join them. 

Let us now consider the last vertex Uk ■ The following two cases are possible: 

(a) there is only one arc that joins Uk with Pk-i(6k) and 

(b) there are two arcs that join Uk with Pk-\{6k)- 

Clearly, the case when Uk is joined with Pk-i{6k) by three or more arcs is 
prohibited because of the absence of cycles in corresponding tree Gk ■ 

In the first case Pk-\(6k) is a connex diagram. One can join this diagram 
with Uk by an arc that have the left foot supported either on grey element or on 
the maximal element of the color group. If there are r colour groups in Pk-i {6k), 
then there are k — r grey elements (see Lemma 3.3). There are r maximal color 
elements. Thus one can choose k elements to put the left foot of the arc. The 
right foot can be put on one of the two elements of Uk- Clearly, |Pfe-i| = dk-i- 
So, the case (a) produces 2kdk-i connected diagrams. On figure 3 we illustrate 
this situation. 




Figure 3: Arc to join Uk with Pk-i- Arrows show possible emplacements of the left foot. 

Not let us pass to the case (b). In this case Pk-i(6k) is splitted in two 
connected diagrams. Let us assume that the left element of Uk is connected by 
an arc with the diagram constructed on j vertices. There are j + 1 possibilities 
to put the left foot of the arc. The right element of Uk joined to the component 
of k — 1 — j elements. There are k — j possibilities to do this. The choice of the 
vertices to produce the component of j elements gives ( fc ^ possibilities. Then 
we get the formula (4.3). It is easy to see that on the way described we obtain 
all the diagarms of T>k- 

Theorem 4.1 is proved. 

Proof of Proposition 4-2. Recurrent relation (4.4) resembles very much the 
recurrent relation determining the moments of the famous semi-circle distribu- 
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tion from random matrix theory 

fc-i 

mfe = y2 rrijirik-i-j ■, m = 1 (4-13) 

0=0 

with a parameter v > [18]. It is also known that the numbers nik — nik(v) 
are proportional to the Catalan numbers Ck ; 

m4( „ ) = „» Cl = (4,4, 

One can easily deduce from (4.13) that 

m k < (2v) 2k . 

Comparing (4.4) with (4.13) taken for v 2 = 2, we conclude that 

h k < 8 fe (4.15) 

for all k > 1. This proves regularity of the generating function h(z) in the 
domain z E Ri/$- 

Now let us show that h(z) verifies equality (4.5). It is easy to see that 
h(x), — l/8<a;<l/8is determined by the differential equation 

"'<*> = T^&r * (0) = 1 - <4 ' 15) 

Indeed, rewriting (4.4) in the form 

fe-i j fe-i 

ftfc = ^2 hj hk-i-j + j:^2 h j hk-i-j, (4.16) 

we multiply both parts of (4.16) by x k and after summation over k > 1, we get 
relation ^ 

h(x)-l = xh 2 (x)+ I h 2 {y)dy. 

Jo 

Passing to the derivatives, we get equation (4.15). 
With the help of the substitution 

xh(x) = ip(x) 

we reduce (4.15) to equation 



x(l - 2-0) 

that gives ipe~ 2 ^ — Cx. Observing that C = 1, we get (4.5). 
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To determine the explicit form of coefficients hk, we use the standard tech- 
nique of the contour integration [17]. First let us note that function ip(z) = zh(z) 
verifies the Polya equation 

iP(z) = ze 2 ^ 

and that the inverse function ip*(w) — we~ 2w is regular in the vicinity of the 
origin. By the Cauchy formula, we have 

/ 1 I^ z h 



2tt i J z k+1 

Changing variables by z — tp*(w), we get dz = (1 — 2w)e~ 2w dw and find that 

1 [l-2w 



^ = V- t ^^e 2wk dw. 

Z7T 1 ./ VJ K 



Then 



.! k k - 2 



(k-l)\ 

and (4.6) follows. Proposition 4.2 is proved. 

Accepting that d\ = 1 and introducing the exponential generating function 



00 ofc-lj 
— 2 «fc k 



^) = E ± iP rfe ' (417) 



fc=i 



and taking into account (4.2), one can write formally for the partition function 
(3.2) that 

hm J-\ogZ n (f3,y)~^D(ge 2 °). (4.18) 

n^oo TIC C I 

The rigorous derivation of this equality will be done in separate publication. 

5 The case of constant edge probability 

Now it is easy to prove the following statement that concerns partition function 
(3.2), where X n is determined by (3.3). 

Proposition 5.1. Given k > I, there exists the limit 

lim \cum k {^X n ) =ff fe V^ fc ;p) = g k V w(S k ;p), (5.1) 

where W and w are polynomials in p involving degrees p k+1 , . . . ,p 2k and p = 
e- 2l3 '/(l - e~ 20 ') with 0' = 0- g. 

Proof. We follow the lines of the proof of Theorem 4.1 because the formulas 
(4.8), (4.9) and (4.10) are still valid in the present case with obvious changes. 
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Regarding (4.10), let us consider the term that corresponds to the trivial parti- 
tion It is easy to see that the formula (4.11) reads as 

W(tt£ } ) = (Ea) k ~Ea k = p k+1 . 
Clearly, relation (4.12) takes the form 

E{iW fci a(* fc ))} ' ' ' V{YT^ k M8 k ))} = P k+1+X > 

and this weight does not vanish as it was before. In particular, another trivial 
partition of the set Uk into k subsets produces the factor p 2k . Relation (4.9) 
completes the proof of the first equality of (5.1). 

The second equality in (5.1) represents the another order of the summation: 
we fix a partition TTk and consider the sum over all diagrams 5k G T>k 

W (5 k ;p) = (-iy-\a-iy. J2 E{f Tl(7rfc;g( , fc)) }---E{f TCT( ^ ;a( , fc)) }. (5.2). 

s k ev k 

Certainly, the numbers W(irk',p) and w(5k',p) are uniquely determined. It is 
possible to obtain recurrent relation that determine the weight W(iTk',p)- This 
recurrent relation generalize (4.3) but it is very cumbersome and complicated. 
We do not present it here. 



6 Discussion and perspectives 

We have introduced the discrete analog of the matrix models with quartic po- 
tentials. We have shown that in this approach the Erdos-Renyi ensemble nat- 
urally arises and plays the same role as that played by GUE for the matrix 
models. Regarding cumulant expansion of the normalized partition function, 
we have shown that the connected diagrams on two-valent vertices replace the 
four-valent two- vertex maps seen in analogous situation in matrix models. Let 
us discuss our results with respect to the properties of random graphs. 
Given a graph 7 with adjacency matrix A(i), the variable 

1 1 " 1 n 

-X n = - £ {A% = - £ An At j (6.1) 

i,j = l i,l,j=l 

represents the number of all possible two-step walks over 7. 

From this point of view, it is natural to ask the same question about asymp- 
totic behavior of the g-step walks and study the terms of the cumulant expansion 
of variable 

n 

e(")(c; 9 ) = -logE^{e^ ) }, with X™ = £ (A") ijt q>2. (6.2) 

It is not hard to show that in the situation of (6.2), the reasonings of Sections 3 
and 4 remains true with obvious changes. The first modification is that instead 
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of the two-valent vertices u s we have to consider q-valent vertices. More pre- 
cisely, u s are represented by the linear graph with q + l subvertices and q edges 
that join them. But it is not hard to see that the number of corresponding dia- 
grams is equivalent to the number of diagrams constructed on the set of g-valent 
vertices, or in other words g-stars with ordered (or labelled) off-spreads. Let us 
count the number of connected reduced acyclic diagrams on q-stars. 

We start with the case q = 3. It is not difficult to repeat the proof of 
Proposition 4.2 and to derive recurrent relations for k > 2 

4 3) = 3(2fc-l) 41 1+ 31 {k > 3} x | ^H-T^'i + + 1 )4?4 3) 



+h^} x ( E zqM® 1 + 1)(2j2 + 1)(2j3 + ^M^h 

Kh+h+h=k-i,ji>i J1 - J2 - Ji - y 

(6.3) 



(3) 

with the initial condition d\ } = 1. Introducing the auxiliary numbers 

hf = {2k + l)df/k\ 

and setting = 1, it is not hard to show that (6.3) is equivalent to the 
following recurrent relation 

e ^s>e *>i- (6-4) 

ji+h+33=k—l,ji>0 



Using (6.4), we derive differential equation 
dh^(x) _ 3(h^(x)) 3 



dx 



l-6x(h( 3 )(x)) 



2 ■ 



h (3 Ho) = i, 



where h^(x) = J2 k >o h k )xk - 

In the general case of q > 2, we obtain equation 

f^M = «(»'""»' r , „««,) = !. (6.5) 

l-(g 2 - g )a;(/i(9)(a;)) 9 

Substitution = x {h^(x)) q 1 leads us to equation 
/(*) = 1 • 



X 1 — ((J 2 — q)ip(x) ' 



Resolving it and returning back to the function h(x) — h^ q \x), we arrive at the 
equation that generalizes the Polya equation 

h(x) = exp{< ? x/i^ 1 (x)}, q>2. (6.6) 
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Using again the function tp and repeating computations of Section 4, we obtain 
relations to determine the coefficients ,k>0; 

E ^-V^fei-tf-gf ^^f ho = l. 

jl+ — +jg-l=k 

In particular, for q = 3, we obtain the first values = 3, = 45/2, 
hf 1 = 1071/6 that correspond to the sequence — 1, — 9, d^ = 153. 

As we have seen, the discrete matrix model (2.5) we proposed represents an 
interesting source of questions about corresponding combinatorial structures. 
The limiting expressions we compute represent the first, zero-order approxima- 
tion to the free energy per site of this model. It could be interesting to further 
develop the diagram approach to study the next terms of the 1 /n-expansion of 
this free energy and justify equality (4.18). A special attention is to be paid for 
the limiting transition n —> oo, c = const. One can expect that the combina- 
torial structure of the cumulants of X n is determined again by the diagrams of 
the type T>j but now all of the diagrams with j = 1, . . . , k are to be involved. 
This is a subject of a separate publication. 
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